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ABSTRACT: 



This program computes the sum of products of deviations of X and Y 
about their means. The computations are performed with double 
precision fixed point arithmetic. 
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1. Title: Double Precision Sum of Squares or Products Subroutine 
Author: Richard A. Lamm 

Date: October 9, 1959 

Installation: Lederle Laboratories, Pearl River, New York 

Classification: F 

2. Purpose: To compute the sum of products of deviations of X and Y 

about their means using the formula: 

Sum of Products = Z XY - (ZX ZY)/n 

Computations are performed with double precision fixed 
point arithmetic. Sums of squares are obtained as the 
special case in which Y = X. 

3. Restrictions: 

a. All data (coded) must be entered at q = 30. The X's 
must be stored sequentially, beginning in Lo(X), and the 
Y's must be stored sequentially starting in Lo(Y). 

b. All data (coded) must be non-negative if a sum of pro- 
ducts is to be generated. If only sums of squares are 
to be generated the data need not be non-negative. 

c. Data coding, when practical, should be such that the 
majority of the data words are four decimal digits in 
length. This recommendation is contingent upon the 
anticipated degree of variation. 

d. ZX, ZY £ 230 _ 1# 

e. In general, the aboslute value of the resulting sum of 
products must be less than 230 - 1 to be of utility. This 
restriction is not applicable if the program using this 
subroutine contains double precision arithmetic. 

f . The double precision word structure is shown on the flow 
chart. H referes to the most significant part of the word, 
L to the least significant part. 

g. No external storage is required for this routine. 

h. All constants and intermediate and final results are 

stored in Lo + 0236 to Lo + 0263. The resultant sum of 
products is designated by Db (DHb and DjLb). 

k. Method: The sum of products is obtained from the standard formula: 

Sum of Products = ZXY -(ZX ZY)/n 

The double precision word is held in two storages such 
that the 29 least significants bits are in one word and 
the remaining most significant bits are in the second 



Title: Double Precision Sum of Squares Program No. - FI-I6I4. 
or Products Subroutine October 9, 19p9 



word. In this way, the 30th position of the most sig- 
nificant part coincides with the first position of the 
least significant part. This Is shown in the word 
structure on the flow chart. 

5. Coding Information: 

a. Storage. Three tracks of memory for the program, constants 

and temporary storage. 

b. Calling Sequence. 

Location Order Address 



«* 


R 




Lo 


«<+ 1 


U 




Lo 


*C + 2 


M 




Lo(X) 
to(Y) 


<X + 3 


M 




oc+ k 


Z 




n at q = 29 


*+ 5 


Z or 


U 


xxxx 


oC+ 6 


etc. 




with Sum of Rrod. in ace. 



o( + 5 should contain an error halt (Zxxxx ) if no provision 
is made for double precision arithmetic or a Uxxxx to the 
double precision arithmetic. An exit too* + 5 indicates 
that the Sum of Products is greater than 23°- 1 . 

c. Input. Data must be pre-stored in Lo(x) et seq. and in Lo(Y) 
et seq. 

d. Output. If exit is by o( + 6 the Sum of Products is in the 
accumulator at q » 30. If exit is by ©(+ 5 DHb and DLb are 
stored in Lo + 026l and Lo + 0262 respectively. 

e. Timing. Approximately 2 + 0.5n seconds. 

f. Urogram Stops. 

Reason Action 

Z X * 23° Recode X's 

EY^ 2 5 ° Recode Y*s 

E XY > 2 59 Recode X's and/or Y*s 
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